Spin Susceptibility in Underdoped Cuprates: Insights from a Stripe-Ordered Crystal 
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We report a detailed study ol the temperature and magnetic-field dependence of the spin sus- 
ceptibility for a single crystal of Lai.87sBao.i25Cu04. From a quantitative analysis, we find that 
the temperature-dependent anisotropy of the susceptibility, observed in both the paramagnetic and 
stripe-ordered phases, directly indicates that localized Cu moments dominate the magnetic response. 
A field-induced spin-flop transition provides further corroboration for the role of local moments. 
Contrary to previous analyses of data from polycrystalline samples, we find that a commonly- 
assumed isotropic and temperature-independent contribution from free carriers, if present, must 
be quite small. Our conclusion is strengthened by extending the quantitative analysis to include 
crystals of La2- :c Ba a ;Cu04 with x = 0.095 and 0.155. On the basis of our results, we present 
a revised interpretation of the temperature and doping dependence of the spin susceptibility in 
La 2 - I (Sr,Ba) ;l: Cu04. 
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I. INTRODUCTION 

A central controversy in the field of high-temperature 
superconductivity concerns the nature of the magnetic 
response in the normal state: does it come from lo- 
cal moments or mobile conduction electrons? The bulk 
spin susceptibility, x S i nas been central to the discus- 
sion, but the analysis so far has been ambiguous. For 
example, in underdoped cuprates one experimentally ob- 
serves a decrease in the normal-state x s on cooling below 
some characteristic temperature T max ji When this be- 
havior was first identified in underdoped YBa2Cu306+. T 
through measurements of the Y Knight shift by Alloul 
and coworkers^ it was interpreted as evidence for the 
development of a pseudogap in the electronic density 
of states. Such a perspective attributes the bulk mag- 
netic response to paramagnetism of the conduction elec- 
trons (Pauli paramagnetism), and it is often cited 
In contrast, Johnston^ proposed that \ s (T) consists of 
two-components: 1) a temperature-dependent portion 
that scales with doping and that is associated with 
antifcrromagnetically-correlated Cu moments, and 2) a 
Pauli contribution that happens to be roughly propor- 
tional to the doped hole density. This alternative ap- 
proach also remains popular ,2iSiS 

Clearly, the temperature dependence of % s is closely 
tied to the issue of the pseudogap^ a continuing conun- 
drum and challenge. The nature of the spin response is 
also part of the ongoing discussion on the dynamic sus- 
ceptibility, where the momentum- and energy-dependent 
structure observed by neutron scattering has generally 
been interpreted cither as the response of conduction 
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or 



electrons scattering across the Fermi surfaco 
as fluctuations of local Cu momentsj 14 i 15 While there is 
much debate on the dynamic response, we have seen lit- 
tle critical discussion of the alternative interpretations of 
the static susceptibility, or reconsideration of the assump- 
tions that went into analyses first developed two decades 
ago, typically applied to measurements on polycrystalline 
samples 

In this paper, we present bulk susceptibility 



measurements performed on a single crystal of 
Lai.875Bao.i25Cu04. We go through a stepwise analysis 
of the data, considering the extent to which we can dis- 
tinguish between contributions from mobile charge carri- 
ers and local moments. One criterion for such a distinc- 
tion is anisotropy: Pauli paramagnetism is expected to 
be isotropic, whereas the paramagnetic response of local 
moments depends on the orientation of the applied field 
due to anisotropic gyromagnetic factors. Thus, from the 
fact that we observe a temperature-dependent anisotropy 
of the susceptibility in the paramagnetic regime, we can 
conclude that there arc significant contributions from lo- 
cal moments. 

Our most direct evidence for local moments comes 
from the temperature regime below the spin-ordering 
transition that has been characterized by neutron 
diffractioni 16 i 17 i 18 There we see an evolution of the 
anisotropic susceptibility that can be simply interpreted 
in terms of the classic expectations for a non-collincar 
antifcrromagnct (see Fig. 1). Furthermore, on following 
the magnetization vs. applied field, we find evidence for 
a spin-flop transition at H = 6 T. Such a transition is 
consistent with an anisotropy of the supcrcxhange inter- 
action between local moments. 

To describe the doping dependence of the suscep- 
tibility of polycrystalline La2- 3; Sr. E Cu04, Johnston* 
proposed that there should be a component x Cu 
due to antiferromagnetically-coupled local moments. 
Furthermore, he assumed that it had the form 
X^ x (x)f(T/T max (x)), where / is a normalized "univer- 
sal" function. He found that he could get a good fit 
to the data if he included a temperature-independent 
component with a magnitude roughly proportional to 
x. The same scaling analysis was later applied to a 
more extensive set of measurements by Nakano et ulZ 
The temperature- independent component is generally as- 
sumed to correspond to the Pauli susceptibility, \ F \ as ~ 
sociated with the doped carriers. 

This successful, but ad hoc, decomposition has been 
widely applied, and we initially attempted to use it on 
our own results to separate out the contribution from 
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FIG. 1: (color online) Schematic of x° u vs T for ( a ) tne 
stripe phase, and (b) a collinear AF order. (- - -) x C " 
in the absence of order, (c) X c U gL/gc, xfioo]dlv/ gib, and 

Xav" = fxpLoo]ffav/Saf, + §X? U Sav/Sc for Lai .875Ba . 125 0uO4. 
TLt marks LTO<->LTT transition, Tso and Too the spin and 
charge stripe order. Lines show ^ 2DHAF f or J = 100 meV, 
including schematic dependence in ordered phase. 

local moments. We found, however, that a x P compo- 
nent of the commonly accepted magnitude is not com- 
patible with the large anisotropy observed in the single 
crystal data. We have strengthened our case by extend- 
ing the analysis to include data for La2- a; Ba.i : Cu04 with 
x = 0.095 and 0.155. We conclude that any temperature- 
independent x P component must be an order of magni- 
tude smaller than that reported in Refs. 

As a result of the much reduced x P > the extracted x 
is larger at low temperature than previously assumed. 
Based on our new insights, we discuss the implications for 
the evolution of magnetic correlations in La2_ 3; Sr :I ;Cu04 
and La2-2;Ba x Cu04 vs. doping and temperature. In par- 
ticular, we show that the maximum of x Cu {T), which 
occurs at T max (:c), is inversely proportional to the effec- 
tive superexchange energy J determined by Raman and 
neutron scattering. 

The rest of this paper is organized as follows. The next 
section describes the experimental methods. In Sec. Ill, 
we present the data and work step-by-step through the 
analysis. The implications are discussed in Sec. IV, and 



a brief conclusion is given in Sec. V. 

II. EXPERIMENTAL METHODS 

The Lai.875Bao.i25Cu04 crystal studied here, with a 
mass of 0.6 g, was cut from the same single- crystal rod 
grown by the travelling-solvent floating-zone method that 
was the source of samples for previous work i 17 i 18 i 19 For 
the complementary measurements on x — 0.095 and 
x = 0.155, similarly sized crystals were used. The 
magnetization M was measured with a SQUID (super- 
conducting quantum interference device) magnetometer 
at Brookhaven (H max = 7 T), and a vibrating sample 
magnetometer at IFW-Dresden (H max = 15 T). In the 
case of Lai.875Bao.i25Cu04, x = M/H was recorded for 
H || [001] (c-axis) as well as H || [100], [010] and [110] (aft- 
plane) in zero-ficld-coolcd (ZFC) and field-cooled (FC) 
modes. [See the high-temperature-tetragonal (HTT) unit 
cell in Fig. Hfb) for definition of directions.] Measure- 
ments of X[ioo] an d X[oio] produce identical results, while 
X[iio] shows clear distinctions. Bulk superconductivity 
(SC) in this sample appears below ~ 2.3 K in H = 20 G. 
SC fluctuations are detectable up to ~ 40 K for H || 
C j 18 i 19 but are suppressed down to ~ 20 K (~ 5 K) by a 
field of 2 T (7 T). 

III. DATA AND ANALYSIS 

Figure [2]^a) shows the raw FC data for \c and X[ioo] 
in Lai.875Bao.i25Cu04. While there is no signature of 
the 235-K structural transition (Tht) from the HTT to 
the low-temperature-orthorhombic (LTO) phase, there 
are clear changes at Tlt = 54.5 K, where the crystal 
transforms into the low-temperature tetragonal (LTT) 
phase. 18 ' 20 We know that charge-stripe order sets in at 
T L t ( I£ i I2ii2i2L22 ^ c increases rapidly below T L t and is 
independent of the magnetic field [except at low T, due 
to diamagnetism from SC fluctuations; see Fig. \3^c) and 
Ref. [l^|. In contrast, X[ioo] only shows a small step at 
Tlt, followed by a significant decrease below 42 K. This 
decrease is not due to diamagnetism from SC, but is 
connected to the onset of spin stripe order at Tso, as 
detected by /^SR22 and neutron diffraction. 16 ' 18 Below 
H ~ 0.5 T, Tso is marked by a clear kink in X[ioo] an d 
Yji mj, a s is shown in the insets of Fig. [3] as well as in 
Ref.ES 

A. Contributions to the measured susceptibility 

The measured susceptibility Xi(T) {i = c, aft) can be 
written as 

Xl (T) = xl(T) + X cmc + xJ V - (1) 

The core diamagnetism ^ corc amounts to —1.01 x 10~ 4 
emu/mol, as determined from standard tables^ The Van 
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Vleck susceptibility of the Cu 2+ ions, \ii depends on 
the direction of the field, but is independent of tempera- 
ture; we will come back to it shortly. We can write the 
spin contribution as 

X-(T) = X? U (T) + X P . (2) 

Here x? u is the contribution from local moments X 

x Cu = N^kj2(ShSi), (3) 

r 

where N is the density of Cu atoms, is the spin com- 
ponent in direction i for the Cu atom at site r; in the 
paramagnetic phase, xf u remains anisotropic, as g a b and 
g c typically differ by 15% for Cu 2+ compounds^ The 
isotropic Pauli susceptibility can be expressed as 2 ^ 

X p = (W2) 2 p(e F ), (4) 

where g 2 and /o(£f) is the electronic density of states 
at the Fermi energy, cf- (In principle, the paramagnetic 
response of conduction electrons due to their spins is par- 
tially offset by their orbital response, in the form of Lan- 
dau diamagnetism.) For a Fermi liquid, the density of 
states is independent of temperature, so that x P is a con- 
stant. As the cuprates do not follow conventional Fermi- 
liquid behavior ) 8 i 10 i 27 x P might be temperature depen- 
dent. 

In simple metals, such as Na or Al, all of the spin re- 
sponse comes from delocalized conduction electrons. The 
doped cuprates are unusual metals, but one can never- 
theless consider the extent to which the charge carriers 
also determine the magnetic response. Thus, we start 
our analysis with the question: Can we explain all of 
the temperature dependence in our single-crystal Xi[T) 
simply in terms of x P (T)? 

To answer this, we consider the anisotropy between 
Xab{T) and x c (T). From Eqs. (P) and ©, we sec that 
there are two sources of anisotropy, xY V an d X? n - Look- 
ing at Fig. [UJa), we see that \c — X[ioo] varies with 
temperature over the entire measurement range, and es- 
pecially below Tso- As noted above, the Van Vleck 




FIG. 2: (color online) (a) X obtained in FC mode with H = 
1 T applied along [100] and [001] in Lai.875Bao.i25Cu04. (b) 
HTT unit cell with definition of directions. 



component is temperature independent, and hence the 
temperature-dependent anisotropy must come from the 
local-moment contribution. Thus, we immediately con- 
clude that the temperature dependence of the spin sus- 
ceptibility cannot come uniquely from a temperature de- 
pendence of the electronic density of states. 



B. Initial analysis of anisotropy 

In principle, we can determine Xi from the anisotropy 
of the measured susceptibility; however, to do this, we 
need to know the values of x^ V an d the gi factors. For 
isolated Cu 2+ ions, one has 

X r = 8/4/eo, (5) 

X V J = 2/4/er, (6) 

g c = 2 + 8A/eo, (7) 

g ab = 2 + 2A/e x , (8) 

where A is the spin-orbit coupling, and the energies eo 
and ei are the crystal-field splittings between the 3d x 2_ y 2 
and the tctragonally-split t 2g orbitals i 25 i 28 In the CuG<2 
planes, the 3d states turn into bands, and one must aver- 
age over the Brillouin zone to properly evaluate these 
quantities. The only calculation^ of such factors for 
cuprates that we know of was performed for hypothetical 
SC2CUO4. While we cannot apply those results directly 
to our case, we can use them as guidance. 

We find that we must resort to fitting the data in order 
to determine the anisotropy factors. To do so, we will 
first consider the case x P = 0; we will later come back 
to consider the impact of a finite x P - From Eqs. J]}-©, 
we find that, in the paramagnetic regime, 

xS72a6 = Xc7<7c, (9) 

or, equivalently, 

(xa b - x coro - x£)lsL = (xc - x coro - xD/gl (10) 

Thus, we must vary the anisotropic factors in order to 
get a unique Xi I dl from the experimentally measured 
Xab and Xc . We note from Eqs. ©-([5]) that 

(g c - 2)/(g ab - 2) = X V 1x2 = 7, (H) 

where one expects 7 < The analysis is relatively in- 
sensitive to the precise value of 7, so we reduce the num- 
ber of adjustable parameters to two by setting 7 = 3, ap- 
proximately equal to the ratio calculated for Sc2Cu04i 2 ^ 
Applying the analysis to the temperature range 100 K < 
T < 300 K, we obtain g ab = 2.19, g c = 2.58, x2 = °- 09x 
10 -4 emu/mol and xY V = 0-28 x 10 -4 emu/mol (see top 
row of Table 1), as well as the polycrystalline averages 
g av = 2.32 and Xa/ = 0.15 x 10" 4 emu/mol. Finally, 
the results obtained for x? u 5av/5c an d xfwo] / 9lb arc 
plotted in Figs.[]Jc) and[3]Jc). 
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FIG. 3: (color online) x of Lai.875Bao.i25Cu04 obtained in FC 
mode, (a) for H || [100]. (b) for H || [110]. Small deviations 
between curves (within error bar at bottom) for different fields 
and magnetometers due to experimental error were corrected 
by shifts in x so that curves match for T > Tlt where no field 
dependence was observed. Insets: FC and ZFC data for H = 
0.1 T showing clear kinks at Tso and Tlt- (c) Comparison of 
the field dependence of X?"~9lv/9i- 
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FIG. 4: (color online) M(H) of Lai.875Bao.i2sCu04 obtained 
in ZFC mode, (a) for H || [100]. (b) for H [| [110]. Small 
experimental errors were corrected so that at 15 T M/H is 
consistent with x m Fig- El ( see error bar). For the 10 K 
curve, the data is limited to H > 3 T to avoid effects of 
superconductivity at lower fields. Inset: derivative dAI/dH, 
calculated after 3-point averaging of M(H) (except for T = 
20 K with H || [100], where 2-point averaging was used due 
to the larger point spacing) . 



C. Spin anisotropy in the ordered state 

Let us consider how the anisotropy is connected with 
the spin ordering. In the classic case of collincar AF 
order, Fig. [lib), onc measures xf| f° r H || S and 
for the two orthogonal orientations; xf u ^ as 
T — > Oji while x5" g rows relative to the paramagnetic 
state. In contrast, for Lai.875Bao.i25Cu0 4 Xc U behaves 
like x^ u j but Xnoo] an d Xmio] eacn seem to correspond to 
l(xf u + x5 u ) ; as indicated in Fig. [T{a). This anisotropy 
is consistent with the stripe model. The stripes rotate by 
7r/2 from one layer to the next, following the crystalline 
anisotropyr22, and the effective spin anisotropy should ro- 
tate as well, with the Cu spins ordering within the ab- 
plane^ (Note that twinning could also yield the same 
average non-collinear response for a completely collincar 
AF.) 



Our fitted values for the Van Vleck susceptibilities are 
quite comparable to the theoretical results^ of Xab = 
0.15 x 1CT 4 emu/mol and x^ V = 0.4 x 10" 4 emu/mol. 
We note that considerably larger values for Van Vleck 
susceptibilities have been reporte d 30 i 31 ; however, those 
studies did not include the possibility of anisotropy in 
xf u - In fact, we will reproduce this effect in Sec. IV. F. 



D. Spin-flop transition 

Next, consider the field dependence of X[ioo] shown in 
Fig. GHa). Increasing the field above 4 T causes Xfioo] 
to rapidly increase for T < 40 K; when H — 15 T, 
X[wo](9lv/glb) matches x? u (5av/5c)- This behavior is 
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FIG. 5: (color online) Model for spin structure of site-centered 
stripes as a function of field. (•, o) Half-filled charge stripes. 
Stripes in adjacent planes at z = and 0.5 are perpendicular, 
(a) spin structure for H = 0, (b) H || [100], (c) H || [110]. 
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FIG. 6: (color online) _ff-T-diagram of La2- a; Ba :I Cu04 as de- 
termined for H || [100]. Only for H < 0.5 T and H > 7 T an 
anomaly at Tso is observed. For SC properties see Ref. OH. 



Lai.8Euo.2Cu04,— suggesting that the in-plane gap has 
changed only moderately. 



consistent with a spin-flop in the spin stripes. The spin- 
flop transition is clearly seen in the Mr 10 o] (H) curves 
of Fig. 0Ja). M[wo] becomes nonlinear for T < 50 K, 
with the transition marked by a peak in dM[ 10 ^/dH at 
Hsf ~ 6 T, as shown in the inset. Similar data for 
X[iio] m Fig-EHb) and M[ 110 ] in Fig. [31(b) exhibit a much 
broader transition with no peak in dM\no]/dH [see inset 
of Fig. [5Ja)]. The anisotropy in dM/dH implies that the 
spins that flop at the transition initially have S || [100]. 
Figure [5] illustrates the situation. At zero field, we pre- 
sume that S || [010] in plane z = 0, and S || [100] in plane 
z = 0.5. For H > Hsf along [100], spins in plane z = 0.5 
flop so that the staggered moments are approximately 
perpendicular to the field [Fig. [5](b)] . It is the gain in 
Zeeman energy, by slightly polarizing the spins in the 
field direction, that stabilizes the flopped state. When 
applying H || [110], spins in all planes continuously ro- 
tate until the staggered moment again is perpendicular 
to H [Fig. EJc)] , so there is no sharp transition. 

Figure [5] shows the resulting phase diagram for H | 
[100]. The LTO^LTT transition at T LT is field in- 
dependent, while the onset of spin stripe order shifts 
from T so = 42(3) K at 0.1 T to 49(3) K at 15 T. X- 
ray diffraction data indicate that charge stripes order at 
Too = T LT M2L21 Hence, for T so < T < T co there 
are static charge stripes but no static spin stripes, which 
is consistent with recent inelastic neutron scattering re- 
sults indicating the opening of a spin wave gap only for 
T < T S o^ On the other hand, \ c increases right at Tco, 
which suggests that spins become parallel to the CuC>2 
planes; see Fig. (21a). This means that below Too spin 
dimensionality is effectively reduced from 2D-H to 2D- 
XY. (Exchange anisotropy becomes relevant even in the 
absence of static order.) This is extremely interesting be- 
cause a 2D-XY system allows for topological order with- 
out additional anisotropics or interlayer coupling, with 
2D SC being one possibility^^^ Finally, Fig. © indi- 
cates the spin-flop transition for T < Tso- -^sf is of the 
same magnitude as in the LTT phase of antifcrromagnctic 



E. Diamagnetism in paramagnetic state 

The reader may already have noted in Fig. [3Jc) the 
apparent anisotropy at T > 54 K, where X? u 9av/9c < 
xfioo] dlv / 9 ab- This behavior does not make sense in 
terms of exchange anisotropy. Instead, we believe that 
X^ u bears a weak two-dimensional diamagnetic contri- 
bution associated with the vortex liquid state proposed 
by Li et al.^ based on magnetization measurements on 
La2- 2; Sr 2; Cu04 for 0.03 < x < 0.07. Similar behavior ap- 
pears above T c at other compositions of La2- x Ba 2; Cu04. 
Examples for x = 0.095 and 0.155 can be seen in Fig. [7] 



F. Anisotropy and doping-dependence of 
susceptibility 

Now we will reconsider our initial assumption that 
\ p = 0. Obviously, that assumption is in conflict 
with the scaling analysis , 6 i 7 i 31 which gives approximately 
X P ~ x. To properly evaluate the latter form, we need 
to jointly analyze anisotropic measurements for samples 
with different hole contents. In Fig. [DJa) we show the 
raw data of Xi on single crystals of La2- 2: Ba :l; Cu04 with 
compositions from x = 0.095 to 0.155. Keeping in mind 
that xY^ and gi should be independent of doping, we 
now repeat our fitting of these parameters for all three 
samples simultaneously and with various choices for x F ■ 
These choices and the fitted parameter values are listed 
in Table I, and the best-fit results for xf n are plotted in 
Fig. [ZUb-d) . The quality of fit parameter R is equal to the 
sum of the squared differences between the two versions 
of Xi (ffav/fff )' with H || ab and H \\ c, evaluated over 
the fit interval of 100 K < T < 300 K. 

Case I corresponds to x F = 0, and clearly it gives an 
excellent match of Xi {dt-v 1 9i ) f° r a ^ three compositions. 
The fitted parameters are quite close to what we obtained 
for x = 0.125 alone (case 0). For cases III and IV, we fix 
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TABLE I: Parameters from analysis of \l f° r different cases, as discussed in the text. Bold numbers indicate parameters 
that were varied during fit. Numbers in the first row result from individual fit for crystal with x = 0.125; see Fig. QIc). In 
contrast, the data in rows I-IV were obtained by simultaneous fits for all x, see Fig. [7Jb-d). Susceptibilities are given in units 
of 10~ 4 emu/mol. R is the sum of squared differences. 



Case 


p 

Xx=0.095 


p 

Xx=0.125 


Xx=0.155 


Xab 


xr 


Sab 


9c 


R 












0.09 


0.28 


2.19 


2.58 


0.005 


I 











0.11 


0.33 


2.18 


2.56 


0.16 


II 


0.3 


0.3 


0.3 


0.18 


0.54 


2.18 


2.56 


0.16 


III 


0.66 


0.89 


1.13 


0.11 


0.33 


2.58 


3.74 


36 


IV 


0.66 


0.89 


1.13 


0.35 


1.04 


2.03 


2.10 


9.4 



X P (x) to values interpolated from results in Refs.HQH^- 
In case III, we use the x^ V °f case I an d only vary the gi, 
while all parameters are varied for case IV. In both, the 
matching is very poor. In fact, for case IV the resulting 
X? u (9lv/9i ) becomes negative as a result of the large 
anisotropy in x7 V that is required to compensate for the 
reduced magnitude of xf n {9%, 1 9i) ■ With reduced values 
of xJ V i as m Case III, the gi factors become unreasonably 
large. 

It is important to note that the analysis is insensi- 
tive to a doping-independent but finite % p . Although 
unphysical, we present this situation in case II, which 
demonstrates a match of equal quality to Case I. We in- 
clude this case simply to show a limitation in terms of 
uniqueness of our analysis. 

While the values of x P (a;) inferred from Refs. 
clearly seem to be unreasonable when anisotropy is taken 
into account, it is still possible that there is a x P { x ) 
contribution of smaller magnitude. Hence, we have per- 
formed fits allowing x F ( x ) ~ x with a variable scale fac- 
tor. Figure[7Je) shows R as a function of the correspond- 
ing x P ( x = 0.155). We see that R strongly increases 
for x F ( x = 0.155) > 0.25 x 10 -4 emu/mol, which is a 
limit nearly five times smaller than the value given by 
the scaling analysis (see case IV in Table I) . We conclude 
that if a temperature-independent contribution from mo- 
bile carriers is present, it must be an order of magnitude 
smaller than the estimates based on the scaling analysis 
of the static magnetic susceptibility. 



IV. DISCUSSION 

Our analysis indicates that the bulk spin susceptibility 
X s in Lai s75Bao.i25Cu04 comes dominantly from local 
Cu moments and that x P does not make a substantial 
contribution. In Fig. 1(c), wc compare our results for x Cu 
with the prediction for the 2DHAF with J = 100 meV, 
where the estimate for J is taken from a neutron scatter- 
ing study of the magnetic excitation spectrumJi A clear 
implication of our analysis is that x s of the doped cuprate 
is greater than that of the undoped antiferromagnet. A 
comparison of the anisotropy below Tso suggests a sub- 
stantial ordered moment, as found by muon spin-rotation 



studiesi 23 i 38 If we estimate x^ n as 



xff u = 2x{ioo](^ v /3^) - x^KJai), (12) 
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FIG. 7: (color online) (a) Xi( T ) of La 2 - I Ba a: Cu04 for dif- 
ferent dopings x — 0.095, 0.125, and 0.155. Note that the 
crystals with x = 0.095 and 0.155 are bulk superconductors 
with T c — 32 K and 30 K, respectively, (b-d) Spin susceptibil- 
ity X? u (T)(gav/9i) extracted from data in (a) under different 
conditions, as discussed in the text. Roman numbers I — IV 
refer to corresponding parameter sets in Table [T] (e) Least 
squared deviation R vs \ P at a; = 0.155 doping for fits in 
which a linear x dependence of x F ( x ) is assumed. Open cir- 
cles indicate R values for cases I and IV from Table I. 
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600 
T (K) 

FIG. 8: (color online) Xav" vs T for the Lai.875Bao.i25Cu04 
single crystal and for polycrystalline La2- 3; Sr a ;Cu04.— Solid 
lines correspond to ^- 2DHAF f or different J.— 



then it appears that x^ u remains finite at low tempera- 
ture, in contrast to the expectation xf| for an un- 
doped antiferromagnet. The lack of perfect spin order- 
ing was already indicated by a Cu nuclear quadrupole 
resonance (NQR) study^ where attempts to measure 
the antifcrromagnctic resonance at temperatures down 
to 0.35 K provided evidence for a broad distribution of 
hypcrfine fields at the Cu site (averaged over the NQR 
time scales). 

For another comparison, Fig. [H] shows Xav" data for 
polycrystalline La2- 2; Sr a; Cu04 with x = 0.08 and 0.15, 
extracted from the bulk susceptibility using the same 
values for the core diamagnctism and Van Vleck sus- 
ceptibility as for the Lai.875Bao.i25Cu04 result, which 
is shown as a polycrystalline average. We note 
that the La2- x Ba x Cu04 result is consistent with the 
L^-zSr^CuC^ data both in magnitude and in terms of 
the trend of increasing Xa V u with doping. We note that 
the latter observation is contrary to what has been in- 
ferred based on the scaling analysis^ 

In Fig. [5J we also compare calculations for the 2- 
dimensional Hcisenbcrg antiferromagnet (2DHAF) with 
different values of J. As one can see, x 2DHAF provides a 
reasonable estimate of Xmax = X Cu (Tmax) for an appro- 
priate J, but gives poor agreement for the temperature 
Tmax at which x Cu reaches its maximum. Now, we know 
from neutron scattering studies that doping significantly 
modifies the dynamic spin correlations from those of the 
undoped insulator ) 17 ' 41 ' 42 ' 43 ' 44 so it should not be sur- 
prising that T max would be reduced from that of a spin- 
only model system. On the other hand, Xmax mi § 

ht be 

less affected, as it corresponds to the response as anti- 
ferromagnetic correlations just begin to overcome distur- 
bances from thermal or electron-scattering excitations. 
Indeed, Imai et al£^ concluded quite some time ago, on 
the basis of NQR measurements of the 63 Cu nuclear spin- 
lattice relaxation rate, that the Cu spin correlations in 
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(color online) Comparison of experimental ' 7 ' 37 ' 40 
) (evaluated with our Xa^X) an d values calculated 
from experimental J values^ 7 -^-^ (with / = (A/2z eB )N = 
0.094iV for z c g =4), as discussed in the text. 



FIG. 9: 

Xav ax 



La 2 _ 2 ;Sr :l; Cu04 are independent of doping at high tem- 
peratures. 

Johnston^ has collected analytical results for a wide 
range of low-dimensional Hcisenbcrg models (with only 
nearest-neighbor coupling) and has shown that Xmax has 
the general form 



Xav (^max) 



.4 



-off 



J 



(13) 



where z c g is the effective number of nearest neighbors 
and A is a scale factor that ranges from 0.4 for a system 
of spin dimers (z e g = 1) to 0.75 for a 2D square lat- 
tice (z e f{ = 4). Thus, if z e & does not change much with 
doping, we might expect Xmax to vary as N e s/J, where 
Aoft = (1 — x)N and N is the density of Cu atoms. 

In Fig. [51 we plot x m ax extracted from the 
literatur o 6 ' 7 ' 37 ' 40 and analyzed with our scheme. For 
comparison, we have used Eq. (fl~3|) . with z e ff = 4 and 
J taken from neutro n 17 ' 41 ' 42 and Raman 4 ^ studies, to 
estimate Xmax- For the case x P = 0, the experimental 
and calculated trends for xSax are quite similar across 
the underdoped regime. Clear deviations appear as x 
crosses into the overdoped regime. Of course, T max — ► 
for x ~ 0.22 f 7 ^ 7 - suggesting that we get to a point 
beyond which dominant antifcrromagnctic correlations 
never fully develop. [However, x(^) remains temperature 
dependent at high doping, and at least up to x ~ 0.3 does 
not demonstrate the temperature-independent behavior 
expected for a conventional metal.—] Our interpretation 
is consistent with neutron-scattering observations of a 
rapid fall off in the weight of antifcrromagnctic spin fluc- 
tuations with overdoping . 48 ' 49 For comparison, we also 
show the case Xmax = Xmax — X P obtained from the scal- 
ing analysis ^2^L which is clearly inconsistent with our 
expectation based on Eq. (CO 
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The developing picture in vmderdoped cuprates is one 
of considerable incoherence of the electronic correlations 
at high temperatures.— Optical conductivity studies of 
I^-zSr^CuO^ I^-zBazCuOzi, and YBa2Cu306+ K in- 
dicate that the Drude peak, often taken as an indica- 
tion of coherent behavior, disappears as the tempera- 
ture is increased towards the pseudogap crossover tem- 
perature, j 1 *, 51 ] 52 ! 53 A recent analysis 5 ^ of temperature- 
and doping-dependent Hall-effect measurements^ 5 - sug- 
gests that not only the doped holes but also thermally- 
activated carriers participate in the transport at high 
temperature. It has been note d 54 ' 56 that T* identified 
from transport properties is quite close to T max from 
the bulk susceptibility. Thus, it appears that the co- 
herent nodal metallic states^ 3 - and the antiferromagnetic 
spin correlations develop together on cooling below T*. 
This co-evolution is compatible with the development of 
stripe correlations j 5 - 7 - as suggested by neutron scattering 
studiesJ&S 8 .^ 

There have been suggestions that the pseudogap 
might be the consequence of spin-density-wave (SDW) 
correlations ! 60 ' 61 In that case, the magnetic response 
would come dominantly from the charge carriers. This 
possibility now seems unlikely, as we have argued that the 
static susceptibility is dominated by the response of local 
moments, and there is no evidence of coherent electronic 
states at T > T* from which a SDW might develop. 

Finally, our interpretation of the static susceptibility is 
compatible with the distinct dynamic responses detected 
by nuclear magnetic resonance studies 6 ^ at Cu and O 
in-plane sites in La2-xSr.j;Cu04. The site-dependent re- 
sponses can be understood in terms of the spatial inho- 
mogeneity associated with stripe correlations! 64 ' 65 ' 66 ' 67 ' 68 
By symmetry, the O site does not see the spins on the 
two nearest-neighbor Cu sites in the uniform antiferro- 
magnetic state; however, stripe correlations can break 



this symmetry, allowing O to detect the spin suscep- 
tibility. In fact, it is the existence of the stripe cor- 
relations that allows for the coexistence of regions of 
antiferromagnetically-coupled Cu moments and mobile 
charge carriers! 18 ' 57 



V. CONCLUSION 

To conclude, we have presented evidence that the 
bulk susceptibility in underdoped and optimally-doped 
La2- 2; (Sr,Ba) a ;Cu04 is dominated by the response of 
antiferromagnetically-coupled Cu moments, with little 
contribution from free carriers. This conclusion has im- 
plications for the interpretation of the dynamic mag- 
netic susceptibility. If electronic quasiparticles were 
to contribute significantly to the dynamic susceptibility 
through Fcrmi-surface-nesting effects, then one would ex- 
pect an associated contribution to the static susceptibil- 
ity. The absence of a substantial isotropic component of 
electronic origin in x s raises questions about the relative 
importance of quasiparticlc contributions to the dynamic 
susceptibility in underdoped cuprates. 
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